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ABSTRACT 

HD 189733 is a K2 dwarf, orbited by a giant planet at 8.8 stellar radii. In order 
to study magnetospheric interactions between the star and the planet, we explore the 
large-scale magnetic field and activity of the host star. 

We collected spectra using the ESPaDOnS and the NARVAL spectropolarimeters, 
installed at the 3.6-m Canada- France- Hawaii telescope and the 2-m Telescope Bernard 
Lyot at Pic du Midi, during two monitoring campaigns (June 2007 and July 2008). 

HD 189733 has a mainly toroidal surface magnetic field, having a strength that 
reaches up to 40 G. The star is differentially rotating, with latitudinal angular velocity 
shear of d£t = 0.146 ± 0.049 radd -1 , corresponding to equatorial and polar periods of 
11.94 ± 0.16 d and 16.53 ± 2.43 d respectively. The study of the stellar activity shows 
that it is modulated mainly by the stellar rotation (rather than by the orbital period 
or the beat period between the stellar rotation and the orbital periods). We report no 
clear evidence of magnetospheric interactions between the star and the planet. 

We also extrapolated the field in the stellar corona and calculated the planetary 
radio emission expected for HD 189733b given the reconstructed field topology. The 
radio flux we predict in the framework of this model is time variable and potentially 
detectable with LOFAR. 

Key words: stars: magnetic fields - stars: planetary systems - stars: activity - stars: 
individual: HD 189733 - techniques: spectropolarimetry 



1 INTRODUCTION 

Magnetic fields are present at different scales in the uni- 
verse, from planets, to stars, galaxies and galaxy clusters. 
Thanks to new high-resolution spectropolarimeters, we are 
able to study the large scale magnetic field of stars, trying 
to understand its origin, but also its implication in different 
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stellar phenomena (stellar wind, magnetic braking, stellar 
cycles, ...). In Hot Jupiter (HJ) systems (giant planets or- 
biting close to their parent stars, i.e. with semi-major axis of 
the planet lower than 0.1 AU), the study of the stellar mag- 
netic field is important to understand Star-Planet interac- 
tions (SPI), and their effects on the evolution and properties 
of the system. 

SPI can be of two types: magnetospheric (e.g. caused by 
reconnections between the stellar and the planetary mag- 
netic fields), or tidal (r esulting from the proximity and 
masses of the two bodies). ICuntz et al.l (|2QQQh suggested that 
such interactions may enhance t he stellar activity. Studying 
the a ctivity of HJ hosting stars, Shkolnik et al . (2003. 120051 . 
2008) concluded that not all the observed systems show hints 
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of interactions and that for a single system with interactions, 
those interactions may be not observable during some ob- 
serving epochs, yet observable at other times. 

Different theoretic al scenarios of magnetospheric inter- 
actions were proposed. IPreusse et al.l (120061) descri bed SPI 
by adopting the Alfven wind model, lLanzal (|2008l ) consid- 
ered a non-potential magnetic field configuration for the 
close d corona of the s t ar. To explain the 'on-off ' nature of 
SPI, ICranmer Saarl (|2007h studied their signatures over 
many orbital cycles considering a cyclic stellar magnetic 
field. These signatures do not repeat exactly from orbit to 
orbit, nor from epoch to epoch. SPI is an intermittent phe- 
nomenon depending strongly on the configuration of the stel- 
lar field. 

The study of the stellar magnetic field is therefore im- 
portant to understand SPI. We have started an observing 
program aimed at detecting and modelling the magnetic 
field of HJ hosting stars. Eleven systems are observed, hav- 
ing different stellar and planetary parameters. In this paper, 
we present the results for HD 189733. This system is inter- 
esting to study SPI, with a short orbital period different 
from the stellar rotation period. 

The properties of HD 189733 are listed in section [2] In 
sections [3] and [U we will present our data, data analysis and 
magnetic modelling of the star. Results of its magnetic topol- 
ogy and differential rotation will also be presented. Stellar 
activity will be studied in section [5] From magnetic maps 
of the stars obtained in section [U we will extrapolate the 
magnetic field in the stellar corona (section [6} , deduce the 
expected planetary radio emission (section [7]) , and end with 
our conclusions (section [8}. 



2 PROPERTIES OF HD 189733 

HD 189733 is a well-known p l anet h osting star, discov- 
ered in 2005 by iBouchv et all (|2005h . The K2V star is 
bright, nearby and active (V=7 1 7 j T eff = 5050 =b 50 K, 
[Fe/H] = -0.03zb0.04 lBouchv et al1l2Q05l d=1 9.3 pc). It has 
a mass of 0.82 ± 0.03 (IBouchv et al.1 2005|i and a radir 
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R* = 0.76 ± 0.01 R (jWinn et alJl2QQ7h . IWinn et al.1 (|2006h 
measured vsini (2.97 ± 0.22 krns" 1 ). The rotation period 
of th e star was measured by different teams using photome- 
try. iHebrard &; Lecavelier Pes Etangsl (|2006l ) found it to be 
about 11.8 d. while IWinn et al.l ( 20071 ) found quasi-periodic 
flux variations of 13.4 d period, which they attributed to 
the stellar rotation. These different values may be due to 
temporal variations in spots coupled to differential rotation 
over stellar latitudes. 

As the planet transits HD 189733, all the planetary 
parameters are well constraine d. A recen t revi ew of the 
planetary orbit was done bv iBoisse et al.l ([2009) . The in- 
clination angle of the orbit is of 85.76° =b 0.29°, the planet 
has a mass of 1.13 ± 0.03 and a radius of 1.154 ± 

0.032 R^ . It orbits the star every 2.2185733 ± 0.0000019 d 
on a circular orbit. T he semi-majo r axis of the orbit is 
of 0.031 ± 0.001 AU ([Boisse et al.1 (|2009h and references 
therein) . iTriaud et al.l ( 20091 ) measured the projected spin- 
orbit misalignment angle and found it to be 0.85° 

— 0.28' 

both orbital and stellar rotation axes are almost aligned. 



3 OBSERVATIONS 

We observed HD 189733 in June 2007 and July 2008 us- 
ing ESPaDOnS and NARVAL spectropolarimeters, the for- 
mer installed at the 3.6-m Canada-France-Hawaii Telescope 
(CFHT) and the latter at the 2-m Telescope Bernard Lyot 
(TBL) at Pic du Midi. These spectropolarimeters provide 
spectra with a resolution of 65000 which span the whole op- 
tical domain (370 to 1000 nm). Each polarization spectrum 
is extracted from a sequence of four subexposures, taken 
in different configurations of the polarimeter waveplates, in 
order to perform a full circular polarization analysis. 

The data were reduced using a fully automatic reduc- 
tion tool, called Libre- Esp rit, installed at CF HT and TBL 
for the use of observers ([Donati et al.l ll997). The spectra 
have error bars at each wavelength pixel and are normalized 
to a unit continuum. Their wavelength scale refers to the he- 
liocentric rest frame. They are automatically corrected from 
spectral shifts resulting from instrumental effects (e.g. me- 
chanical flexures, temperature or pressure variations) using 
telluric lines as a reference. Though not perfect, this pro- 
cedure allows spectra to be secured w ith a radial velocity 
(RV) precision of better than 30 ms -1 (jMoutou et alJl2007l : 
iMorin et alJl2008h . 

We collected 20 spectra in June/ July 2007 using both 
ESPaDOnS and NARVAL. These data covers about 2 stellar 
rotations (25 nights); NARVAL data cover roughly the first 
rotation cycle, and ESPaDOnS data sample well the second 
rotational cycle. In 2008, we only used NARVAL. Our data 
are spread over 14 nights, from 10 to 24 July. This second 
data set covers only slightly more than one rotation cycle. At 
around 700 nm, the S/N ratio of the spectra varies between 
220 and 940 for 2007 spectra, and from 390 to 720 for 2008 
spectra (per 2.6 k ms -1 velocity bin). W e also use the 2006 
data published in iMoutou etaD J200 7) merged in a single 
data set (as opposed to lMoutou et al.l (|2007h who analyzed 
separately the June 2006 and August 2006 data) . The phase 
coverage is poor for this epoch. The observations, separated 
by 52 days, overlap in terms of rotation cycle. The complete 
log of t he 2007 and 2008 o bservations is given in Table [T] 
(refer to lMoutou et al.ll2007l . for information about the 2006 
data) . 

The rotational and orbital phases, denoted ER Q t and 
Eorb, were computed using the two ephemerides: 

To = HJD 2, 453, 629.389 + 2.218574 E 0rb 

To = HJD 2, 453, 629.389 + 12 E Rot (1) 

The first ephemeris is that of IMoutou et al.1 {2007). For the 
ephemeris giving the rotation phase, we use a rotation pe- 
riod of 12 d, identified as the equatorial rotation period (see 
section [4} . 

To improve the S/N ratio of our data and extract the 
polarization from many lines simultaneously, we used Least- 
Squares Deconvolution (LSD) which assumes that all lines 
more or less repeat the same polarization information. In 
practice, it consists of deconvolving the observed spectra by 
a line mask, computed using Kurucz's lists of atomic line 
parameters (Kurucz CD-Rom 18) and a Kurucz model at- 
mosphere with solar abundances, temperature and logarith- 
mic gravity (in cms -2 ) set to 5000 K and 4.0 respectively. 
The line mask includes the moderate to strong lines present 
in the optical domain (those featuring central depths larger 
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Table 1. Journal of June 2007 and July 2008 observations. Columns 1-11 sequentially list the UT date, the instrument used, the 
heliocentric Julian date and UT time (both at mid-exposure), the complete exposure time, the peak signal to noise ratio (per 2.6 kms -1 
velocity bin) of each observation (around 700 nm), the rotational and orbital cycles (using the ephemeris given by Eq. [T}, the radial 
velocity (RV) associated with each exposure, the rms noise level (relative to the unpolarized continuum level I c and per 1.8 kms -1 
velocity bin) in the circular polarization profile produced by Least-Squares Deconvolution (LSD) and the longitudinal magnetic field. 



Date (UT) 


Instrument 


HJD 


UT 




S/N 


Rot. Cycle 


Orb. Cycle 


^rad 


^LSD 


Bi 










(2,454,000+) 


(h:m:s) 


(s) 




(52+) 


(284+) 


(kms" 1 ) 




(G) 


09 


June 


07 


NARVAL 


260.532660 


0:42:06 


4x900 


210 


0.5953 


0.4817 


-2 


213 


1.59 


-3.5 + 3.2 


13 


June 


07 


NARVAL 


264.561950 


1:23:52 


4x900 


390 


0.9311 


2.2978 


-2 


382 


0.84 


-0.2 + 1.7 


18 


June 


07 


NARVAL 


269.559040 


1:19:11 


4x900 


330 


1.3475 


4.5502 


-2 


152 


1.01 


-1.0 + 2.0 


19 




07 


AT A T>T T A T 

NARVAL 


271.479030 


OO OO A O 

23:23:48 


a v , c\r\r\ 

4x900 


360 


1.5075 


5.4156 


-2 


299 


0.97 


— 1.3 ± 1.9 


23 


June 


07 


ESPaDOnS 


275.004850 


12:00:39 


4x900 


920 


1.8013 


7.0049 


-2 


294 


0.34 


-2.1 ± 0.7 


23 




07 


ESPaDOnS 


275.121020 


14:47:56 


4x800 


870 


1.8110 


7.0572 


-2 


340 


0.36 


-0.6 ± 0.7 


26 


June 


07 


ESPaDOnS 


277.847340 


8:13:36 


4x900 


740 


2.0382 


8.2861 


-2 


416 


0.45 


-3.2 ± 0.9 


26 


June 


07 


ESPaDOnS 


278.117560 


14:42:43 


4x900 


870 


2.0607 


8.4079 


-2 


352 


0.36 


-3.3 ± 0.7 


27 


June 


07 


ESPaDOnS 


278.841160 


8:04:38 


4x800 


910 


2.1210 


8.7341 


-1 


997 


0.36 


-6.2 ± 0.7 


28 




07 


ESPaDOnS 


279.840490 


8:03:35 


4x900 


810 


2.2043 


9.1845 


-2 


390 


0.40 


-2.9 ± 0.8 


30 


June 


07 


ESPaDOnS 


281.863400 


8:36:26 


600+140 


270 


2.3729 


10.0963 


-2 


325 


1.19 


-2.0 ± 2.3 


01 


Tulv 


07 


ESPaDOnS 


282.834040 


7:54:05 


4x700 


770 


2.4538 


10.5338 


-2 


164 


0.44 


-2.2 ± 0.9 


01 


Julv 


07 


ESPaDOnS 


283.120060 


14:45:56 


4x800 


690 


2.4776 


10.6627 


-2 


039 


0.49 


-3.7 ± 1.0 


02 


Tulv 


07 


ESPaDOnS 


283.949800 


10:40:42 


4x600 


590 


2.5467 


11.0367 


-2 


286 


0.57 


-5.1 ± 1.1 


02 


Tulv 


07 


ESPaDOnS 


284.131990 


15:03:02 


4x600 


640 


2.5619 


11.1188 


-2 


331 


0.52 


-3.9 ± 1.0 


03 


Tulv 

j uiy 


07 


ESPaDOnS 


284.975470 


11:17:36 


4x600 


760 


2.6322 


11.4990 


-2 


262 


0.42 


—3.5 ± 0.8 


03 


Tulv 


07 


ESPaDOnS 


285.112190 


14:34:28 


4x600 


660 


2.6436 


11.5607 


-2 


149 


0.50 


-2.9 ± 1.0 


04 


Tulv 

j uiy 


07 


ESPaDOnS 


285.842200 


8:05:38 


4x600 


510 


2.7044 


11.8897 


-2 


089 


0.65 


-5.0 ± 1.3 


04 


Tulv 


07 


ESPaDOnS 


286.014210 


12:13:19 


4x600 


490 


2.7188 


11.9672 


-2 


167 


0.68 


-3.6 ± 1.3 


04 


Tulv 

j uiy 


07 


ESPaDOnS 


286.087050 


13:58:13 


4x900 


450 


2.7248 


12.0001 


-2 


184 


0.74 


-3.9 ± 1.5 


10 


July 


08 


NARVAL 


658.455050 


22:47:44 


4x900 


570 


33.7555 


179.8412 


-1 


977 


0.55 


-3.7+1.1 


11 


July 


08 


NARVAL 


658.603660 


2:21:43 


4x900 


590 


33.7679 


179.9082 


_2 


014 


0.56 


-2.7+1.1 


11 


July 


08 


NARVAL 


658.636700 


3:0:9:18 


2x900 


360 


33.7706 


179.9231 


_2 


025 


0.90 


-4.4 + 1.8 


14 


July 


08 


NARVAL 


662.448330 


22:37:52 


4x900 


440 


34.0883 


181.6411 


-1 


970 


0.74 


-3.8 + 1.5 


15 


July 


08 


NARVAL 


662.589210 


2:00:44 


4x900 


480 


34.1000 


181.7046 


-1 


942 


0.69 


-1.4 + 1.4 


15 


July 


08 


NARVAL 


662.633270 


3:04:11 


4x900 


510 


34.1037 


181.7245 


-1 


939 


0.65 


-1.8 + 1.3 


15 


July 


08 


NARVAL 


663.446250 


22:34:51 


4x900 


360 


34.1714 


182.0909 


-2 


230 


0.95 


1.8 + 1.9 


16 


July 


08 


NARVAL 


663.587830 


1:58:42 


4x900 


440 


34.1832 


182.1548 


-2 


284 


0.78 


3.3 + 1.5 


16 


July 


08 


NARVAL 


663.631880 


3:02:09 


4x900 


390 


34.1869 


182.1746 


-2 


298 


0.90 


4.2 + 1.8 


17 


July 


08 


NARVAL 


665.486420 


23:32:37 


4x900 


730 


34.3415 


183.0105 


-2 


129 


0.42 


-0.3 + 0.8 


18 


July 


08 


NARVAL 


665.587170 


1:57:42 


4x900 


720 


34.3498 


183.0559 


-2 


168 


0.43 


-1.1 + 0.9 


18 


July 


08 


NARVAL 


666.466700 


23:04:11 


4x900 


660 


34.4231 


183.4524 


-2 


193 


0.48 


0.0 + 0.9 


19 


July 


08 


NARVAL 


666.569380 


1:32:03 


4x900 


630 


34.4317 


183.4987 


-2 


122 


0.49 


-0.3 + 1.0 


19 


July 


08 


NARVAL 


666.657390 


3:38:47 


2x900 


420 


34.4390 


183.5383 


-2 


063 


0.70 


1.2 + 1.4 


20 


July 


08 


NARVAL 


668.491560 


23:39:55 


4x900 


690 


34.5919 


184.3651 


-2 


263 


0.47 


-8.8 + 0.9 


21 


July 


08 


NARVAL 


668.590150 


2:01:53 


4x900 


680 


34.6001 


184.4095 


-2 


224 


0.46 


-7.3 + 0.9 


21 


July 


08 


NARVAL 


669.465410 


23:02:15 


4x900 


430 


34.6730 


184.8040 


-1 


939 


0.78 


-5.0 + 1.6 


22 


July 


08 


NARVAL 


669.565650 


1:26:35 


4x900 


510 


34.6814 


184.8492 


-1 


952 


0.63 


-1.5 + 1.3 


22 


July 


08 


NARVAL 


669.659660 


3:41:57 


2x900 


420 


34.6892 


184.8916 


-1 


968 


0.72 


-2.4 + 1.4 


22 


July 


08 


NARVAL 


670.469810 


23:08:33 


4x900 


580 


34.7567 


185.2567 


-2 


337 


0.55 


-2.2 + 1.1 


23 


July 


08 


NARVAL 


670.570050 


1:32:54 


4x900 


670 


34.7651 


185.3019 


-2 


323 


0.47 


0.9 + 0.9 


23 


July 


08 


NARVAL 


670.651980 


3:30:52 


2x800 


430 


34.7719 


185.3389 


-2 


305 


0.69 


1.1 ± 1.4 


24 


July 


08 


NARVAL 


672.478300 


23:20:44 


4x900 


680 


34.9241 


186.1621 


-2 


266 


0.47 


-4.2 + 0.9 


25 


July 


08 


NARVAL 


672.577800 


1:44:01 


4x900 


720 


34.9324 


186.2069 


-2 


277 


0.44 


-4.9 + 0.9 



than 40% of the local continuum, before any macroturbulent 
or rotational broadening, about 4,000 lines throughout the 
whole spectral range) but excludes the strongest, broadest 
features, such as Balmer lines, whose Zeeman signature is 
strongly smeared out compared to those of narrow lines. In 
addition to the intensity and polarization profiles, LSD pro- 
duces a null profile (labelled N) that should contain no polar- 
ization; this helps to confirm that the detected polarization 
is real and not due to sp urious instrumental or reduction 
effects ([Donati et al . 1997). The multiplex gain provided by 



LSD in V and N spectra is of the order of 30 with regard 
to a single line with average magnetic sensitivity, implying 
noise levels as low as 35 parts per million (ppm). A typical 
set of LSD profiles is shown in Fig. [Q for 20 July 2008. 

The radial velocity (RV) of the star is obtained by fit- 
ting each Stokes I profile with a Gaussian profile. The RV 
values we obtain are listed in Tableland are in good agree- 
ment with the expect ations, when using the orbital solution 
of iBoisse et al.l (|2009h once the data are globally shifted by 
-0.06 kms" 1 and -0.15 kms" 1 for June 2007 and July 
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Velocity (km/s) 

Figure 1. Typical LSD profiles shown for 20 July 2008. 
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Figure 2. Radial velocities of HD 189733 derived from June 

2007 (NARVAL: yellow dots and ESPaDOnS: magenta dots) and 
July 2008 (black triangles) spectra as a function of orbital phase, 
with their error bars (20 — 30 ms _1 ) shifted by 0.06 kms -1 and 
0.15 kms -1 respectively (see text for more details). The radial 
velocity model pl otted here (blue fu ll line) corresponds to the 
orbital solution of iBoisse et al.l (|2009T ). 

2008 respectively (see_Fig. [2}. A similar offset was reported 
bv iMoutou et alj (|20Q7h . These variable RV shifts may be 
due to stellar variability, to the presence of a stellar com- 
panion (Bakos et al. 2006), or to a yet unknown cause. The 
wavelength calibration being done relative to the telluric 
lines (giving a precision of about 30 ms -1 ), we believe it 
is very unlikely that these RV shifts could be due to any 
instrumental effect. 

For the magnetic analysis in the rest of the paper, the 
spectra used are the LSD ones, corrected for the orbital mo- 
tion of the planet. 



4 MAGNETIC IMAGING 

4.1 Model description and imaging method 

To reconstruct the magnetic m aps of HD 189733, w e use the 
Zeeman Doppler Imaging code ([Donati et ali l997) in its lat- 



est version. ZDI is a tomographic imaging technique, it in- 
verts series of circular polarization Stokes V profiles into the 
parent magnetic topology, i.e. the distribution of magnetic 
fluxes and orientations. Since the problem is ill-posed, ZDI 
uses the principles of maximum entropy to retrieve the sim- 
plest image compatible with the data. The latest version of 
the code describes the field by its radial poloidal, non-r adial 
poloidal and toroidal components (|Chandrasekharlll96ll ) , all 
expressed in terms of spherical harmonics expansions. This 
has the advantage that both si mple and com plex topologies 
can be reconstructed reliably (|Donatill200"lh . Moreover, the 
energy of the poloidal and toroidal components, or of ax- 
isymmetric and non-axisymmetric modes can be estimated 
directly from the coefficients of the spherical harmonics ex- 
pansion. Given the small value of vsini, the resolution at 
the surface of the star is limited; we therefore truncate the 
spherical harmonics expansion to the five lowest terms only, 
i.e., I < 5. 

Reconstruction proceeds by iteratively comparing the 
reconstructed profiles to the observed one, until they match 
within the error bars (i.e. reduced chi-square Xr ~ !)• Prac- 
tically, the star is divided into 9000 grid cells of similar 
area. The contribution of each grid cell to the reconstructed 
Stokes profiles is calculated, given the RV of the cell, the 
field strength and orientation in it, the location and pro- 
jected area. Summing the contributions of all the grid cells 
yields the synthetic profile. 

For each grid cell, the local unpolarized Stokes I is mod- 
elled by a Gaussian with a full width at half-maximum 
(FWHM) of 7 kms -1 , central rest wavelength of 500 nm 
and effective Lande factor of 1.25. The Stokes V profile is 
calculated assuming the weak-field approximation, i.e. the V 
profile is proportional to the line-of-sight projected compo- 
nent of the field (called longitudinal field and denoted Bi), 
as well as the first derivative of the local I profile. The incli- 
nation angle of the rotation axis of the star with respect to 
the line-of-sight is ~ 85° (see section [2}. 

For a differentially rotating star, magnetic regions lo- 
cated at different latitu des have different angular velocities. 
As in [Fares et al.l ((2009), we consider that the rotation at the 
surface of the star follows 0,(6) = O e q — dOsm 2 (6), where 
0(0) and O eq are respectively the angular velocities at a 
latitude and at the equator, and dO is the difference in 
rotation rate between the pole and the equator. The po- 
sition of the grid cells are calculated by this law relative 
to the median observing epoch. Measuring the recurrence 
rate of the signatures from magnetic regions located at var- 
ious latitudes gives access to the amount of surface shear. 
In practice, for each pair of (f2 eq , dO) in acceptable range of 
values, we reconstruct a magnetic image at a given informa- 
tion content (constant magnetic energy) from the observed 
profiles and get the Xr °f the reconstruction procedure. Fit- 
ting a paraboloid to the Xr values we obtain by this process 
gives the optimum differential rotation (DR) parameters of 
the star. 



4.2 Results 

4-2.1 Differential rotation 

Including DR for the reconstruction of June 2007 map is nec- 
essary to fit the observed Stokes profiles almost down to the 
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0.5 0.52 0.54 



(rad/d) 

Figure 3. Variations of Xr as a function of f2 eq and dQ , derived 
from the modelling of the Stokes V data set for June 2007. The 
outer colour contour corresponds to a 3.4% increase in the x? > 
and traces a 3 a interval for both parameters taken as a pair. 

noise level. We applied our method to measure the DR for 
this data set. The Xr values as a function of Q eq and dQ form 
a well denned paraboloid shown in Fig. [3] The optimal 
DR parameters, for which Xr is minimal, are respectively 
Q eq = 0.526±0.007 radd -1 and dQ = 0.146=b0.049 radd -1 . 
The equator of HD 189733 rotates in 11.94±0.16 d, while its 
pole has a slower rotation of 16.53 ±2.43 d. This corresponds 
to a time for the equator to lap the pole by one complete 
rotation cycle of 43.06 ± 14.44 d (lap time). 

For July 2008, we are not able to measure DR, the 
Xr map featuring no well defined minimum (over the range 
of parameters we explored). Our data have a lower S/N ra- 
tio than those of June 2007, they cover slightly more than 
an equatorial rotation period; they do not contain all the 
information needed for measuring the recurrence rate of the 
signatures of high latitude magnetic features. For this epoch, 
the DR parameters we include in our reconstruction proce- 
dure are those of June 2007 (for the same Xr value), sup- 
posing that the DR did not change between the two epochs. 

For the 2006 data, we are not able to measure the DR 
given the very small phase interval covered both in June and 
August, the resulting % 2 map yielding again a very chaotic 
surface. Nevertheless, adding the differential parameters of 
June 2007 in our reconstruction gives a better fit to the data 
(Xt ~ 1-26 as opposed to 3.9 for solid-body rotation with a 
rotation period of 12 d). 

4-2.2 Magnetic maps 

For June 2007, the reconstructed Stokes V profiles includ- 
ing the DR values we obtained are shown in Fig. [U the data 
being fitted for a Xr value slightly larger than 1. HD 189733 
has an average surface magnetic field of 22 G. The field has 
a predominant toroidal component, contributing 57% of the 
total magnetic energy. The poloidal component of the field 
is mainly non-axisymmetric, 67% of its energy being in 
modes with m > 1/2. The quadrupole, octupole and higher 
orders contribute almost equally to the poloidal energy. Or- 
ders with / > 3 contains 30% of the poloidal energy, while 
70% is in the lower orders. 

In July 2008, the field has a greater average value of 



Table 2. Average magnetic field on the surface of the star, per- 
centage of the toroidal energy relative to the total one, percentage 
of the energy contained in the axisymmetric modes of the poloidal 
component for each epoch of observation. 



36 G. Its strength reaches up to 40 G in some magnetic re- 
gion (see Fig. [5}- The toroidal component contributes 77% 
of the total energy. While the poloidal field is still mainly 
non-axisymmetric, its dipolar component has a stronger con- 
tribution for this epoch. Orders with / > 3 still have a con- 
tribution of about 32% to the poloidal energy. 

When merging June and August 2006 in one map, the 
poloidal field contributes 67% to the total energy. The av- 
erage value of the field is of 33 G. The properties of the 
magnetic field are listed in Table [2] 

For the three epochs, we do not observe a global change 
in the magnetic polarity. The radial field shows positive mag- 
netic regions covering the pole and a concentration of the 
magnetic regions around the equator for all epochs. A drop 
in the poloidal energy is observed between 2006 and 2007, 
as well as between 2007 and 2008 (by a factor of 1.6 and 
1.9 respectively). The topology of the field also changes over 
two years, with the toroidal component strengthening sig- 
nificantly at the expense of the poloidal component. 



5 ACTIVITY INDICATORS 

The Ca n H & K and Ha lines are tracers of the stellar chro- 
mospheric activity. We studied the variability of the resid- 
ual emission in those activity proxies for our data. For each 
tracer, we calculated a mean profile per run (see Fig. [6] in 
the particular case of July 2008) , then subtract it from each 
spectrum. The residual emission profiles obtained for the 
Ca II H &K are shown in Fig. Fitting those profiles with 
a Gaussian gives the residual emission value (the equivalent 
width of the Gaussian, in kms -1 ). 

We studied the variation of the residual emission for 
each observing epoch (including 2006). The star exhibits 
variability on both short (hours) and long (days) time scales. 
We modelled the longer term variability as the sum of a sine 
wave (to describe the periodic modulation) and a linear func- 
tion (to describe the longest term variability) . For a range of 
period between 2 d and the length of the run, the period for 
which the Xr °f the fit is minimal is the best fit period to the 
data. The best fit periods are listed in Table [3] We plotted 
the Ca II H & K residual emission and the best fit solution 
in Fig. [8j For all three epochs, the Ca II H & K residual 
emission is modulated on a time-scale of the equatorial ro- 
tation period. For Ha, the residual emission is modulated 
on a time-scale of ~ 13 d (on average), again close to the 
rotational period of the equator. Ha and Ca n H & K resid- 
uals are well correlated (see Fig. [9} except at 2 main epochs 
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Figure 4. Maximum-entropy fits (thin red line) to the observed Stokes V LSD profiles (thick black line) of HD 189733 for 2007 June 
(left) and 2008 July (right). The rotational cycle of each observation (as listed in Table [T} and la error bars are also shown next to each 
profile. 



Table 3. The best fit period for the Ca II H & K and the Ha 
residuals in days (using the model described in the text). The 
error bars correspond to a 3 a error bars on our measurements. 



Epoch 


Ca II Residual 
Period(d) 


Ha Residual 
Period(d) 


2006 


12-I5±g;l 


n-6t ; 3 2 


June 2007 


12.1+2- 2 


14.4±i;| 


July 2008 


11 9 +1 - 2 


I3.8t!;§ 



We then searched for periodic modulation of these resid- 
uals for a range of periods by fitting the residuals with a 
sine wave. For June 2007, we find a period of 3.65+q'7 5 d 
in Ha, significant to 3 a. We do not find a similar period 
in the Ca n H & K residuals. For July 2008, two periods 
of 2.45+0*25 d and 3.5+o - 2 d were found in Ha and a large 
peak around 2.9 d in the Ca n H & K (2.9±° 5 6 d). Only 
in July 2008 for the Ha emission, one of the two identified 
periods is roughly compatible with the beat period. 



(2006 August 8-12, 2007 June 9-19) where Ra is slightly 
stronger at a given Ca n H & K residual emission. 

To look for additional low-amplitude periodic fluctu- 
ations, we subtracted the rotational modulation from the 
data (the best fit solution obtained). In the case of SPI, the 
emission variability would be modulated by the orbital pe- 
riod or rather by the beat period between the stellar rotation 
and the planetary orbital period (varying between 2.5 and 
2.7 d, depending on the stellar latitude considered). Fig. [TOl 
shows these residuals as a function of the orbital phase for 
the Ca II H & K . A large dispersion in the residuals can 
be seen. In particular, we can see that the residual emis- 
sion is highly dispersed at any given orbital phase range, 
directly reflecting the short term variability of the activity 
proxies. We therefore suspect that this is mainly caused by 
intrinsic variations in the activity of the star rather than to 
SPI, as we would expect enhanc ed emission variabilit y to be 
concentrated at specific phases (|Shkolnik et al.ll2008h . 



6 EXTRAPOLATION OF THE MAGNETIC 
FIELD 

From the surface magnetic field, one can then extrapolate 
the magneti c field in the stellar atmosphere assuming a po- 
tential field (|Jardine et al.l l2002). In the case of SPI, extrap- 
olating the magnetic field lines is a way to understand the 
environment in which the close-in planet evolves, and to in- 
vestigate potential interactions (e.g. reconnection events due 
to interactions of the stellar and planetary magnetic fields). 

The extrapolation technique was first applied to the so- 
lar corona, and then to other stars. As for the Sun, we as- 
sume that there is a surface beyond which the field becomes 
purely radial, named the Source Surface (hereafter SS). All 
closed field lines, connecting regions on the stellar surface 
of different magnetic polarity, are inside the SS (the heights 
of the magnetic loops are smaller than the SS radius). From 
the extrapolation of the magnetic field, we can calculate the 
magnetic energy at any given point inside and on the SS. To 
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Figure 5. Maximum-entropy reconstructions of the large-scale magnetic topology of HD 189733 as derived from our 2006 (top row), 
2007 June (middle row) and 2008 July (bottom row) data sets. The radial, azimuthal and meridional components of the field (with 
magnetic flux values labelled in G). The star is shown in flattened polar projection down to latitudes of —30°, with the equator depicted 
as a bold circle and parallels as dashed circles. Radial ticks around each plot indicate rotational phases of observations. 
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Figure 6. Mean Ca II H & K (left panel) and Ha: (right panel) profiles for July 2008. 
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Figure 7. Call H Sz K residual profiles for June 2007 (left panel) and July 2008 (right panel). The rotational cycle of the star is 
mentioned next to each profile. 




Figure 8. Equivalent width variations in the emission cores of the Ca II H & K (after subtracting a mean profile calculated for the run) 
as a function of the Heliocentric Julian Date for June 2007 (left panel) and July 2008 (right panel). The red curves are the best fit to 
the data (see text for more details). 



calculate the magnetic energy outside the SS, we consider 
that the magnetic flux is constant on radial shells. 

In the system HD 189733, the massive planet is at 
8.8 R*. The position of the SS is not well known. We first 
consider that the SS is located at 3.4 R* (as a lower limit 
given that the SS radius Rss ~ 3 R@ for the much less active 
Sun), i.e. inside the planetary orbit. Figure ITT1 represents an 
extrapolation of the stellar magnetic field within the source 



surface for our two epochs of observations. It shows in par- 
ticular that the magnetic configuration in the stellar corona 
is complex. We find that the magnetic field at the distance 
of the planet has a strength that can reach values up to 
40 mG (i.e. 4000 nT, see Fig. [T2|). The magnetic energy is 
null when the planet crosses in front of the intersection of 
the neutral line with the equator. 

We considered a second case for which the SS is at 
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Figure 9. Ha and Ca II H & K residual emission for the three observing epochs (blue, red and blacks dots are respectively for 2006, 
2007 and 2008 data). 
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Figure 10. Residuals of the Ca II H & K equivalent width after subtracting the rotational modulation, as a function of the planet orbital 
phase for June 2007 (left panel) and July 2008 (right panel). 



5.8 R*, also within the planetary orbit. The mean magnetic 
field at 5.8 i?* averaged over the longitudes drops by a factor 
of 2 relative to the previous case. 



Finally, we consider cases where the SS is larger than 
8.8 R*. When the SS is at 10 R*, we notice that the magnetic 
field at the planet orbit is smaller than that for a SS at 3.4 R* 
by a factor slightly larger than 3 for the mean value (see Fig. 
[T2]) . For a SS beyond 10 R*, the magnetic field at the planet 
orbit does not change dramatically relative to the previous 
case. In all these cases, the value of the field is not constant 
along the planet orbit. 



We thus expect average field values of 4 — 23 mG at the 
distance of the planet. 



7 THE EXPECTED RADIO EMISSION 

Knowing the stellar magnetic field, one can estimate the 

expected radio emission for HD 18 9733b. 

In the magnetic energy model (Zark a et al. I l200lh . plan- 
etary radio emission is believed to be energized by the Poynt- 
ing flux transported by the stellar wind. The amount of 
power emitted by radio waves P ra d is roughly proportional to 
the power input Pmp ut supplied b y the stellar wind magnetic 
energy (see Figure 6 IZarkdl2007h . Thus, the radio power is 
given by 

Prad OC Pinput OC Ve^B^R^ (2) 

where v e ff and B± are respectively the velocity of the stellar 
wind and the component of the interplanetary magnetic field 
(IMF) perpendicular to the stellar wind flow, both in the 
reference frame of the planet, and R s denotes the radius 
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Figure 11. The extrapolated magnetic field of HD 189733 for June 2007 (left) and July 2008 (right). White lines corresponds to the 
closed magnetic lines, blue ones to the open field lines (reaching the source surface). The red line corresponds to the field line joining the 
planet to the star at a given (arbitrarily selected) orbital and rotation phase. 




longitude longitude 



Figure 12. The magnetic field at the distance of the planet as a function of stellar longitude for June 2007 (left) and July 2008 (right) 
for different positions of the SS (the black, red and blue lines are respectively for a SS placed at 3.4, 5.8 and 10 R*). 



of the planetary magnetosphere (known as magnetospheric 
standoff distance). 

To estimate the radio power, one has to know B±, v e ff 
and R s - The y depend on the stellar magnetic field, rotation 
and age (see GrieB meier et al.ll2007l for more details) . When 
a star is older, its wind strength is weaker. Via the stellar 
wind strength (velocity and density) , the stellar age governs 
the size of the planetary magnetosphere, i.e. the cross section 
on which the planet can intercept the energy flux. 

Given the rotation period of 12 d, the stellar age 
is of rsj 1.6 Gy r (bas ed on the formalism presented in 



GrieBmeier et al 



GrieBmeier et al 



Using the stellar wind model of 
, we calculate R s and v e ff, that are 
equal respectively to 3.4 i? p i an et and 332 kms -1 (both con- 



tance of the planet is taken from the results of the extrap- 
olation (see section [6] and Fig. fT2|) . In the particular case 
where the SS is beyond the planetary orbit, our calculation 
includes the three components of the field. 

In this context, we find average radio fluxes of about 7 to 
220 mJy depending on the location of the SS, for a frequency 
range of 0-6 MHz. Values up to 20 MHz are compatible 
within the model uncertainties. We also find this radio flux 
to be time variable (see Fig. 113)) PI 



stant on the planetary orbit). The magnetic field at the dis- 



1 I GrieBmeier et a~ ] (120071 ) have considered different parameters 
for the system, in particular a stellar age of ~ 5.2 Gyr, a rotation 
period of 28 d, and a radial magnetic field of about 2 G at the 
surface of the star yielding a different radio flux. 
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Figure 13. Expected planetary radio flux at Earth as a function 
of the rotational phase of the substellar point, for three different 
positions of the source surface. Top panel: Radio flux for the mag- 
netic configuration of June 2007. Bottom panel: Radio flux for the 
magnetic configuration of July 2008. 



8 DISCUSSION AND CONCLUSIONS 

In this paper, we present a detailed spectropolarimetric 
study of the star HD 189733, host of a transiting giant 
planet. The star was observed at two epochs (June 2007 
and July 2008). Using Zeeman Doppler Imaging, we recon- 
structed the magnetic maps of the star. With a strength up 
to 40 G, the magnetic field is dominated by the toroidal com- 
ponent at both epochs. This component contributes 57% and 
77% to the total energy respectively and is mainly axisym- 
metric. In contrast, the poloidal component is mainly non- 
axisymmetric. Its contribution to the total energy drops 
from 2006 to 2007 and 2008. We will continue monitoring 
this system to study the magnetic evolution on time-scales 
longer than 2 years and look for a potential magnetic cycle. 

HD 189733 rotates differentially and has a latitudinal 
angular rotation shear of dQ — 0.146 =b 0.049 radd -1 ; the 
star has an equatorial period of 11.94 ±0.16 d and a po- 
lar period of 16.53 ± 2.43 d. These values of the equatorial 
and polar periods bracket all published photometric periods 
within the error bars. 

The star is an active star, variable on small time-scales. 
We analyzed the activity residuals in the Ca n H&K and 
Ha. These residuals are periodic and modulated with the 
rotational period of the star. Active regions apparently con- 
centrate around the equator, given the modulation period 
of 12-13 d. A rotational modulation was also found by 
Moutou et all I2007L iBoisse et all 120091 . and IShkolnik et~al] 
2008. We looked for lower amplitude periodic fluctuations, 
the periods we found are different than the orbital and beat 
period (2.5-2.7 d) except in one occurrence (Ha, July 2008) 
where one of the detected roughly matches the beat period. 

To enlarge our knowledge of this system, we studied the 
magnetic field in the stellar atmosphere using the extrapo- 
lation technique applied to the reconstructed surface mag- 



netic maps. We find that HD 189733 has a complex magnetic 
topology for both epochs. Depending on where the SS (the 
surface beyond which the field is purely radial) is, we find 
that the magnetic field at the distance of the planet (8.8 R*) 
is variable throughout the orbit, of the order of 4 — 23 mG in 
average. 

We also estimated the radio flux ex pected from S PI as - 
suming the magnetic scenario model of IZarka etaD (|200lh 
and find it to be of the order of 10 — 220 mJy on average. We 
also predict it to be variable with time on a time scale equal 
to th e beat period (contrary to previous published predic- 
tions IGriefimeier et al.1 [2007) . The radio flux we predict at 



0-6 MHz is potentially det ectable with LOFAR (see fig. 1- 
3 in IGriefimeier et al.l l2QQ7^ in the coming years. Published 
radio observations only report upper limi ts on the planetary 
flux at higher frequencies (307-34 7 MHz. lSmith et al.ll2009l : 
iLecavelier Pes Etangs et aUl2009L 244 MHz and 614 MHz) - 
providing no constraint on the model discussed in the 
present paper. The variability of the planetary radio flux 
with the subplanetary stellar phase will make the distinc- 
tion between the planetary and stellar radio flux in the ob- 
servational data more challenging; however, this effect can 
be used to distinguish between the magnetic energy model 
and other models of interactions. Our result confirms that 
a single observation of a star-planet system is not sufficient. 
Rather, it is important to have multiple observations densely 
sampling stellar rotation, planetary orbit and beat periods. 

The study of SPI is an ongoing effort. Monitoring stars 
at different epochs and through multi- wavelength campaigns 
will help us identify the nature of SPI and the origin of 
their apparent on-off behavior. Studying stellar magnetic 
cycles and comparing results for HJ hosting stars with dif- 
ferent stellar and planetary parameters will enlarge our un- 
derstanding of SPI, as well as stellar magnetism and activity 
in general. 
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